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Abstract. Adaptive-network models are typically studied using deterministic differential equations which 
approximately describe their dynamics. In simulations, however, the discrete nature of the network gives 
rise to intrinsic noise which can radically alter the system's behaviour. In this article we develop a method 
to predict the effects of stochasticity in adaptive networks by making use of a pair-based proxy model. The 
technique is developed in the context of an epidemiological model of a disease spreading over an adaptive 
network of infectious contact. Our analysis reveals that in this model the structure of the network exhibits 
stochastic oscillations in response to fluctuations in the disease dynamic. 
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1. Introduction 

The importance of network structure to the dynamics of interacting agents in the real world is now almost 
universally recognised. Traditionally, systems of interacting agents were modelled by means of simple 
ordinary differential equations, describing the deterministic 'mean-field' behaviour of well-mixed populations. 
Whilst this approach has been very successful in the past (examples can be found in ecology, evolution and 
game theory [3 [5]), it is valid only in simplified models in which each agent is equally likely to interact with 
any other agent. This caveat places a severe constraint on the applicability of these techniques to real- world 
systems. 

In recent years it has become possible to go beyond models of well-mixed populations and to move towards 
the incorporation of network structure. This step forward has been facilitated by several recent advances: 
the increasing quality and availability of data on naturally occurring networks; the development of a more 
comprehensive theory of complex networks (see e.g. [3]); and lastly, the exponential growth in computing 
power allowing for large-scale simulations of network models. In these models, individual agents are placed 
on the nodes of a networlQ and interactions occur across the edges of the network. Broadly speaking there 
are two classes of individual-based models on networks (i) those in which the underlying network is assumed 
to be static and (ii) models which take into account changes to the structure of the network concurrently with 
the agent dynamics. The former class focusses on the dynamics on the network and is a good proxy when 
the evolution of the network structure is slow compared to the dynamics of the nodes. The latter approach 
applies more generally as it considers dynamics both on and of the network, occurring simultaneously. 
This second approach is frequently referred to as the 'adaptive network' modelling paradigm, and it is an 
important recent development, providing a useful framework to model many real- world systems [5]. The 
conceptual simplicity of these models is, however, countered by the fact that it is notoriously difficult 
to develop a succinct mathematical theory for their macroscopic behaviour. The equilibrium properties 
of adaptive- network models have been studied in depth [5], frequently by means of the so-called pair 
approximation (PA), usually attributed to [6]. This approach is widely used in adaptive networks; although 
the motivation and quality of the approximation are questionable [7], results (for equilibria at least) are 
typically reasonable. 



f We are here only concerned with models in which there is one single individual at each node of the network. An example of 
so-called meta-population models, with multiple individuals at each node, can be found in [4]. 



The mathematical analysis of the actual dynamics of adaptive-network models and their approach to 
equilibrium is intricate; in particular, the effects of what is referred to as 'intrinsic noise' or 'demographic 
stochasticity' are not at all well-understood. Such stochasticity can have a profound impact on the dynamics 
of individual-based models; a realisation which carries with it important implications for fields such as ecology 
[5JH] and epidemiology [ID] . A robust theoretical framework with which to describe effects of intrinsic noise 
has been developed for well-mixed models, primarily based around the system-size expansion method of van 
Kampen [llj . but no comprehensive approach exists for models on adaptive networks. 

It is the purpose of this article to propose a general theoretical method to address effects of intrinsic noise in 
adaptive-network models. Starting from the rules of the adaptive-network model, we derive a low-dimensional 
Markov jump process which (approximately) captures the aggregate dynamics of the agents and edges in the 
network model. We refer to this model as the 'pair-based proxy' (PBP). In deriving the PBP we employ a 
moment-closure assumption equivalent to the PA. Our implementation is somewhat different to usual in that 
we are not attempting to derive a set of ordinary differential equations; instead we define a new stochastic 
process designed to reproduce the macroscopic behaviour of the original networked model. In doing so we 
gain instant access to established methods applicable to jump processes. In the limit of infinitely large 
networks the jumps in the PBP become infinitesimally small, and the dynamical rules reduce to a system of 
differential equations. Stochastic effects in large but finite networks can be investigated using the linear-noise 
approximation [llj . 

One particularly important application of adaptive-network models is to the study of infectious disease, 
and we have chosen to develop our method in this context. Epidemiological models involving well- mixed 
populations, or even static networks, are at best a crude approximation to the dynamics of real-world 
epidemics. Adaptive-network models represent a step towards greater realism, capturing the dynamics of a 
disease which spreads over a constantly evolving network of infectious contact between individuals. Moreover, 
adaptive networks offer the possibility to model the effects of intervention strategies aimed at disrupting the 
network of infection, such as contact tracing and quarantining. It is also becoming clear that demographic 
noise has an important role to play in the dynamics of epidemics, particularly the noise-induced excitation 
of transient oscillations, leading to cycles of disease outbreaks [T2l UOj . 

The first adaptive network epidemic model is due to Gross et al. [13] , who studied the equilibrium properties 
of the model by means of the PA (they provide a description of the full phase diagram, including an active non- 
equilibrium phase, together with an analysis of the structural properties of the network). This deterministic 
approach is successful in describing stationary states and phase behaviour, but it systematically neglects 
effects of noise. This work was followed by Shaw and Schwartz [T3] who added extrinsic noise 'by hand' to 
the deterministic equations. The method of Rozhnova and Nunes |15j addresses intrinsic fluctuations due 
to demographic stochasticity, however, it is only applicable to static networks in which each node has the 
same number of neighbours. Our approach is similar to that of [15] . who reverse-engineered a stochastic 
pair-based model to fit the macroscopic PA equations obtained from the network model, but crucially we 
here address the case of adaptive networks. We expect that the general scheme we introduce will be useful 
in quantifying stochastic effects in a wide range of adaptive-network models, not only in epidemiology but 
in other areas as well. As a key result of our analysis we are able to identify oscillations of the underlying 
network structure itself, induced by the stochasticity of the dynamics. The (approximate) mapping onto the 
PBP allows us to make analytical predictions for the spectral properties of these oscillations of the network, 
as well as for those of the population of agents on the network. 

The remainder of this paper is organised as follows: In the next section we define an adaptive-network model 
of an epidemic and show how to derive the low-dimensional PBP model which approximates its dynamics. 
The next section contains a theoretical analysis of the PBP for large networks, in which expressions for 
deterministic dynamics and stochastic corrections are computed. In section 4, we compare the theoretical 
results from the PBP with simulations of the adaptive-network model and discuss some of the limitations 
of the approximation. In section 5 we report on the observation that the structure of the network itself 
undergoes noise-driven oscillations, and show how this behaviour can be captured by the PBP. 
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2. Constructing the pair-based proxy 

2.1. Adaptive network model 

We consider a susceptible-infected-recovered-susceptible (SIRS) disease model in a cohort of N individuals, 
which are placed at the nodes of a network and joined by a total of K edges. The edges of the network 
represent the potential for infectious contact over which the disease may spread: whenever a susceptible (S) 
individual shares an edge with an infected (I) individual, the susceptible individual becomes infected with 
rate ft. Infected individuals become recovered (R) with rate 7, and recovered individuals become susceptible 
again with rate 8. If we label individuals by i = 1, . . . , N, and if € {0, 1} is the adjacency matrix of the 
network (we consider only undirected networks, i.e. = ctji), then the dynamics can be described by the 
reactions 

Si + Ij — % Ii + Ij , Ii — ^> Ri , Ri — > Si, (1) 

where the notation Si indicates that the individual at node i is in the susceptible state (and similar for 7$ 
and R4). 

In the real world, the transmission of a disease requires spatial proximity, and thus the network of potential 
infection is in a constant state of flux caused by the movement of individuals. We incorporate this fact into 
our model through a process of constant low-level re- arrangement. Each edge in the network decays (i.e. is 
removed) with rate fj,, at which point a new edge is generated between a random pair of individuals elsewhere 
in the network^. 

As well as this random rewiring process, the structure of the network adapts intelligently in response to the 
progression of the epidemic by a process of 'smart' rewiring. For each susceptible-infected pair there is a 
chance that the susceptible individual will discover the infection and choose to remove that edge and replace 
it with an edge to a randomly chosen susceptible individual. This happens with rate w for each SI pair. 
When referring to the state of the networked model, we write [A] for the number of nodes in state A 
and [XY] for the total number of edges in the system between pairs of nodes in states X and Y, for any 
X,Y £ {S,I,R}. Later we will discuss the behaviour of the model in the limit of large network size, by 
which we mean that N — > 00 and K — > 00 with N/K held constant. 

2.2. Pair-based proxy model 

As discussed in the introduction, we intend to study the macroscopic behaviour of the adaptive-network model 
through a low-dimension model which captures the main features of the network model whilst remaining 
analytically tractable. We will define a 9-dimensional Markov jump process with state vector 

X = (lSj [/J [i?] ISSj ISIj ISR} fRRj [IRj [I/flf . 

The components of the state vector are intended to represent the corresponding quantities in the network 
model: our goal is to derive the transition probabilities for the jump process from the network model in such 
a way that the statistics of, for example, [5] and in the PBP model will reflect those of [S] and [SI] in 
the original network model. 

For each of the five processes taking place in the network model (i.e. infection, recovery, rebirth, random 
rewiring and smart rewiring) we will compute the average net change to each [X] and [Ay], and use this 
to specify the jump direction and magnitude in the PBP. The jump rates are also drawn from the network 
model in this way. To illustrate this construction, we discuss in detail the process of infection. Our argument 
proceeds along the following steps: 

1. The rate with which an infection event takes place in the network model is /3[SI]. The rate for infection 
events in the PBP model is therefore set to /3[S7]]. 

§ The total number of edges in the network remains constant under this rewiring process. It is also possible to consider edge 
creation and deletion as separate processes so that the total number fluctuates, however, this choice does not make a qualitative 
difference to the aggregate behaviour of the model. 
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When an infection takes place in the network, one 5 individual is converted into an /. For the PBP 
we write this as a jump with direction S I and magnitude 1, meaning that [5] i-> [5] — 1 and 

w^w + i- ' ' " 1 

An infection event in the network model changes the configuration of the system on the level of edges 
as well, and several edges will be affected, namely those between the newly infected individual and its 
neighbours. So an infection event will generally result in one (or several) SI edges becoming II edges, 
and similarly some SR edges may become IR, and some SS may become SI. It is not obvious at first 
sight how many of each of the different types of edges are changed during an infection event. 

For the purposes of the PBP model we use the typical change to [XY] which takes place in the network 
model during an infection event. In order to compute these it is first necessary to determine the average 
number of neighbours an 5-individual has in the network: 

2 [SS] + [SI] + [SR] 



[S] 



(2) 



In the PBP model this quantity is approximated by [fc s ] = (2 [55] + [57] + [5i?])/[5]. 

5. Next we determine how many neighbours of a newly-infected individual we expect to be in states 5, 
I and R respectively. In the network model, the probability that a randomly chosen neighbour of an 
5-individual is of type 5 (or I or R) is 2[SS]/([k s }[S}) (or [SI\/([k s )[S\) or [SR]/([k s ][S\), respectively). 

6. In the network model the states of the individuals surrounding the individual of type 5 are correlated. 
In order to formulate the PBP we assume that only correlations via the central node are relevant. 
This is equivalent to the usual pair moment closure assumption used in the literature [6j 1131 1141 115) . 
Making this assumption, the magnitude of the jump in the direction 55 i— > SI is 2z[55]/[5], where 
z = ([Ms] — l)/[fe]- Similarly, the jumps in directions SI H> II and SR h->- IR have magnitudes 
1 + z [51]/ [5] and z[5i?]/[5], respectively. 

7. Each of the jumps computed above arc taken to happen simultaneously; the combined jump is calculated 
by vector addition. To summarize, infection events in the PBP model occur with rate /3[57] and lead 
to the changes to the state vector X as indicated in the first column of the stoichiometric matrix in Eq. 
© below. 

A full list of the transition events, rates and state changes is shown in Table [T] The column labelled 'jump 
rate' indicates how many events of the different types (infection, recovery, rebirth, smart rewiring and random 
rewiring) occur per unit time in the PBP. For clarity, in the 'jump direction' and 'jump magnitude' columns 
we have split the jumps into individual contributions of the form X H ► Y or XZ H > YZ. Note that the jump 
magnitudes for the edges are typically not integer. 

The effects of infection, recovery, rebirth and smart rewiring are also summarized in the following 
stoichiometric!]] matrix: 

/ -1 1 



T (1) 



-1 



(3) 



1 -I 



i -iswm 

[SI]/[5] 

-2[//]/[I] 
(2p7] - p7?])/pl 

The first column for example indicates that the number of individuals of type 5 is reduced by one in the 
event of an infection, the number of I increases by one, the number of R remains unchanged, and the number 
of particles of type 55 is reduced by 2z[55]/[5], and so on. The second, third and fourth columns of the 
above matrix describe recovery, rebirth and smart rewiring events; the corresponding entries can read off 

|| We are slightly abusing terminology here, as some entries of our 'stoichiometric' matrix are not integers, and moreover they 
depend on the state of the system. This is not the normal convention. 



-2z[55]/[5] 
z(2[55]-[5/])/[5] 
-z[SR}/[Sl 
z[5/]/[5] + l 
z[5i?]/[5] 




ism/m 
im/m 
(2[RR}-ism)/m 



-UR]/m 

-21RR]/IR} 






1 

-1 







/ 
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Event 


Jump rate 


Jump direction 


Jump magnitude 


Infection 




5^7 
SS^SI 
SI i y II 
SR^IR 


1 

2zissyisj 

H-z[SJ]/[S] 
zlSRj/lSj 


Recovery 


7 PI 


I^R 
57^ 577 
II i y IR 
IR^ RR 


1 

urn /in 


xieDirtn 




R^S 
SR^SS 
IR^SI 
RR^SR 


i 

isRi/im 


Smart rewiring 




SI^ ss 


i 


Random rewiring 
{A,B,C,De {S,I,R}) 


H(2 - 5 c .D)lABl{C\iDl/N 2 


AB i y CD 


i 



Table 1. List of reactions in the pair-based microscopic model. Here z = ([fes] — l)/[fes]i where 
[fc s ] = (2[SS] + [57] + [SK])/[S] is the average degree of a susceptible node. 



from the table. The rates with which infection, recovery, rebirth and smart rewiring occur are summarized 
in the rate vector 

" ~~ S{R} ' [ > 

V Hsil J 

The random rewiring process is composed of 36 separate reactions, and the associated stoichiometric matrix 
T^ 2 ) (of dimensions 9 x 36) and rate vector R^ are thus rather large. They are constructed as follows: 
for each combination A,B,C,D e {S, I, R} the corresponding column of T^ 2 ^ has a —1 in position for 
JAB], a +1 in the position for [CTJ], and zeros elsewhere; the corresponding entry of the rate vector is 
yu(2 — < l >c,D)[AB]][C][7)]]/iV 2 . The stoichiometric matrix and rate vector for the whole model are found by 
concatenating TW with T^ 2 ) and R (1 1 with . We write 



T = ( T« TT< 2 ) ) 



R 



rw 



(5) 



where T is now a matrix of dimensions 9 x 40. 

The initial condition for the PBP is taken to match that of the network model (i.e. \S\ = [S], [7] = [I], 
and so on). Note that the sum of the first three entries in each column of T is zero, as is the sum 
of the remaining entries in each column. This fact implies the conditions [5] + [7] + [77] = TV and 
[55] + [57J + [577] + [77J + P"#l + = K (which are inherited from the initial condition of the 

network model) must hold throughout the dynamics. 



3. Analysis of the pair-based model 

3.1. Effective langevin equation 

We now proceed to derive an analytical description of the PBP model in the limit of large TV. The first 
step is to introduce the density vector x = X/N and rescaled rate vector r = R/N. There is no need to 
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adjust the stoichiometric matrix though, as its entries are invariant under scaling of X. The scaled variable 
a; is a realisation of a density dependent Markov process satisfying the conditions of Kurtz' theorem [16] . 
Employing that result, we can immediately write down a Langevin equation for x (to be interpreted in the 
Ito sense) which holds when N is large: 

± x = Tr + 'r,(t), (6) 



where rj is a vector of multiplicative Gaussian noise variables with state-dependent correlator 

40 
k=l 

The variable k here runs through all 40 reactions of the PBP model. This formulation is equivalent to that 
obtained by the performing a Kramers-Moyal expansion [11] of the Master equation and neglecting terms 
including and beyond 1/N. 

3.2. Deterministic limit 

In the limit of large population size (that is, N 1, K 1 with constant K/N) Eq. ([6]) reduces to 
a deterministic system of differential equations given by x = Tr. Returning to the original scaling, the 
deterministic equations are: 

dm 



dt 

dm 

dt 

dim 

dt 



= SIR} - PISIl 
PISIJ - 7 [J], 

■■-rin -sim, 



dt = o\SR\ - 2z(3 — + wlSIj - ulSSj + A/x-p-, 

d ^ J l i SS \ ,\ xm i fen , , UQnx9jr PIP] 



dt V ' IS] (Sj J 1 » " 11 y ™ J * N 2 

« = 7[ 5i] + * mm - phi) - ^™ - ,ism + 2^™, 

^= OT fl + ,@V(27 + ,)[I/] + ^ J12 



*3 = 7 [/H] - 2i[JiJi] - ARR] + K^-. (8) 

Here every term is of order N and lower order terms (including Gaussian stochastic corrections) have been 
neglected. 

This system is similar to that studied in |14j . though there are several differences. Firstly, the random 
rewiring process has not been considered previously, and we also take a slightly different prescription for 
the smart rewiring process. More importantly, in most work using the PA, it is implicitly assumed that the 
average degree of susceptible nodes is simply K/N, not \ks\ as we use here. 

With nine dimensions and seven parameters, the bifurcation structure of our deterministic model is likely to 
be complex. Indeed, analysis of a closely related model [14] reveals bifurcations of the transcritical, saddle- 
node and Hopf types. We are most interested in the presence of a Hopf bifurcation, whereby a stable fixed 
point becomes unstable to create a limit cycle. We can be sure that our model possesses such a transition, 
since in the limit of large /i it reduces to the usual well-mixed SIRS model, which does indeed have a Hopf 
bifurcation. For the remainder of the paper will restrict our attention to regions of parameter-space in which 
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the deterministic model has a globally attractive endemic (meaning that [7] 7^ 0) fixed point x*. We aim to 
provide a description of the effects of demographic noise in this regime. 



3.3. Linear noise approximation 

We next obtain a first-order approximation for the effects of noise following the standard procedure of 
linearising the Langevin equation ([5]) around the endemic fixed point x*. In population-based models this 
procedure is equivalent to the van Kampen system size expansion [11] , The central limit theorem suggests 
we should expect fluctuations of order 1 / vN, motiving a change of variables 



£ = VN(x - x*) . (9) 
We apply this change to Eq. © and retain only leading-order terms in 1/y/N to obtain 

j t Z = A£ + C(*) , (10) 
where A is the Jacobian matrix of the deterministic part of (|6]), that is, 

(11) 



d 



and £ is a vector of Gaussian noise variables whose correlation matrix B has entries 



k 



(12) 



where < Ci(t)Cj >= Bij6(t — t'). If the fixed point is stable and the Jacobian A has complex eigenvalues, 
then the trajectories described by Eq. (|TU)) may exhibit stochastic oscillations. To characterise this behaviour 
we take a Fourier transform, writing 

£>) = -L / me-^dt. (13) 



/2tt 

In frequency space the object of interest is the power spectral density matrix, defined by P(w) = 
f (£(cj')£(u;)' i ') duj', where (•••) denotes the ensemble average. For linear noise processes we have the 
explicit formula 

P(w) = (A-iwI)- 1 B(A T + iwI)- 1 , (14) 

which can be derived directly from Eq. (fTUf and using the properties of the noise £■ The quantity I is here 
the 9x9 identity matrix. The diagonal entries of P(uj) contain the power spectrum of oscillations for each 
of the entries of £. For example, fluctuations in the number of susceptible agents are described by Pi,i(cj). 
The off-diagonal entries of P(w) give information on the correlations between state variables. For example, 
oscillations in the density of susceptible agents are typically out of phase with those of infective agents, 
leading to complex values of Pi^{u) with argument close to 7r/2 [T7] . 



4. Test against simulations 

As with all approximation schemes of this type, the motivation for creating the PBP was to provide theoretical 
predictions that are a reasonable match to the behaviour of the adaptive-network model in simulations. In 
Fig. Q] we show an example of the comparison between theoretical results of the PBP and data gathered 
from a corresponding simulation of the adaptive-network model (performed using Gillespie's algorithm (18)). 
As shown in Fig. [TJ the agreement is generally good, though we should point out that there is a particular 
area of parameter space in which there are significant discrepancies between the models. When both the 
mean degree K/N and the rate of random rewiring \i are small, the typical states of the adaptive-network 
model are highly heterogeneous across the network. This heterogeneity makes a description of the network 
only in terms of pairs too weak to capture the details of the model, and the agreement with the PBP breaks 
down, particularly in the transient part of the dynamic. 
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Figure 1. Left: time series of the fraction of susceptible agents coming from solution of the PBP 
deterministic equations J8J (thick dashed blue line) and a simulation of the adaptive network (thin solid 
black line). Right: Power spectrum of oscillations in the fraction of susceptible agents as computed for the 
PBP, see Eq. il 1 41) (thick dashed blue line), and measured from simulations of the adaptive network (thin 
solid black line). The parameters used were f3 = 6, <S = 3, 7 = 0.5, fi = 5 w = 30, N = 10000, K = 100000. 
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Figure 2. Discrepancy between the PBP and adaptive network simulations, as a function of random 
rewiring rate fi. The left hand plot displays the relative error between the equilibrium fraction of susceptible 
individuals in the PBP and adaptive network simulations. In the right hand plot a measure for the relative 
error in the power spectrum of oscillations of susceptible individuals is computed by integrating (from u) = 
to u = 10) the difference between the theoretical result of the PBP and the spectrum measured from the 
network simulations, and dividing by the total power of the simulation result in that frequency domain. The 
other model parameters were (arbitrarily) chosen as = 6, S = 3, 7 = 0.4, ui = 4, TV = 1000, K = 10000. 



On a mathematical level, the root cause of the difficulty in approximating adaptive-network models with 
low-dimensional systems is the development of long-range correlations within the network. Heuristically, it 
is easy to imagine that the random rewiring process we have introduced will act to curtail the development 
of long-range correlations and move the behaviour of the model closer to that of a well-mixed population. 
Whilst the derivation of the PBP did not involve anything so exact as a formal expansion in we expect 
that the PBP and the adaptive network will agree in the limit fi — > 00. We have tested this prediction by 
varying /1 over two orders of magnitude and measuring the agreement between network simulations and the 
PBP theory; the results are shown in Fig. [2j 

5. Oscillations of the network structure 

It is well-known in the literature on epidemic modelling that demographic stochasticity can excite periodic 
fluctuations in the prevalence of the disease, leading to cyclical outbreaks (see, for example, [4l H0l [T2l 115] ). 
It is natural then to ask if the structure of the network of infectious contact also undergoes stochastic 
oscillations in response to the stochasticity of the disease dynamic. 
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Figure 3. Left: time series of the average degree of susceptible agents coming from solution of the PBP 
deterministic equations J8J (thick dashed blue line) and a simulation of the adaptive network (thin solid 
black line). Right: Power spectrum of oscillations in the average degree of susceptible agents as computed 
for the PBP (thick dashed blue line), and measured from simulations of the adaptive network (thin solid 
black line). The parameters used were the same as those in Fig[T] namely, = 6, 6 = 3, 7 = 0.5, fi = 5 
w = 30, N = 10000, K = 100000. 



This information is contained in the oscillations of the mean degree of susceptible nodes, [ks], in the adaptive- 
network model, which have amplitude of the order of 1/yN. The power spectrum for [fcs] (the PBP model's 
equivalent of [ks]), can be derived from that of the other variables of the PBP. Recalling that 

2[SSj + [Si] + {SR} 2x± + x 5 + x 6 



ks IS] 

we linearise around the deterministic fixed point value kg to obtain 

2£ 4 + £ 5 + g 6 - kfc 



(k s - k* s J 



Carrying out the Fourier transform is straightforward, and gives 

2& (w) + &(w) + &M - fcKiM 



ks(lo) 



•TV 



(15) 



(16) 



(17) 



and thus the power spectrum of oscillations in ks is computed as 

(|£sM| 2 ) = (^*) |(fc5) 2 -Pl,lH+4 J P4,4H+P 5 ,5H+-P 6 ,6M 

+ 4Re[P 4 , 5 (o;)] + 4Re[P 4 , 6 M] + 2Re[P 5fi (uj)] 
- k* s URje[P 1A (uj)] +2Re[Pi l5 (w)] + 2Rc[P lfi (uj)] 



(18) 



In Fig. |3] we show a comparison between the deterministic dynamics and stochastic oscillations in [fcg] in 
the PBP to those of [ks] in simulations of the adaptive-network model. The parameter values are the same 
as those used in Fig [1] Whilst the agreement in both plots is not quite as close as in the comparison of [5J 
to [S] , the PBP still clearly provides a strong indication of the oscillations in network structure. 
It is interesting to explore how the expression of stochastic oscillations in network structure varies with 
parameters. The colourmaps in Fig. [4] show the position and height of the peak in the spectral density, as 
the infection rate ft and smart rewiring rate w are varied over two orders of magnitude. Several interesting 
features are visible in the figures. At low values of ft there is a region (marked A in the figure) in which the 
disease typically dies out and thus no fluctuations are present. Two behaviours are possible in the active 
phase. In region B, stochastic oscillations are present, but the spectrum is dominated by a peak at the 
origin which represents the excitation of low frequency modes (although other, smaller, peaks elsewhere in 
the spectrum are still possible). Outside of regions A and B the power spectrum has its global maximum at 
a non-zero frequency; here the oscillations generally become faster as ft is increased and more powerful as 
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Figure 4. Colourmaps showing the position and height (on a logarithmic scale) of the peak in the power 
spectrum of oscillations in [fcg], as a function of infection rate and smart rewiring rate w, both varying 
on a logarithmic scale from 0.1 to 10. The other parameters are <5 = 3, 7 = 0.2, /1 = 2, K/N = 10. In region 
A the trivial (i.e. disease extinction) steady state is attractive and no fluctuations are present. Region B 
denotes the area of parameter space in which the stochastic oscillations are present, but the spectrum is 
dominated by the zero mode; time series from this region are difficult to distinguish from pure white noise. 



w is increased. In principle, the transition in the power spectrum from being dominated by a peak at the 
origin to a peak at non-zero frequency uj could happen in one of two ways: either the peak at zero could 
migrate, or a secondary peak could grow elsewhere in the spectrum and eventually overtake the one at the 
origin. Since the boundary of region B is sharp, we can conclude that in our model the second mechanism 
is at work. 



6. Conclusions 

To summarize, we have put forward an approximative modelling approach with which to study the effects of 
intrinsic fluctuations in interacting agent systems on adaptive networks. Our method follows a similar line 
to that of [15] : the key step consists of deriving an 'effective' stochastic process for the average numbers of 
single and pairs of agents in each possible state. This technique was previously seen to be viable for static 
regular networks, we here show that it can be successful for the case of co-evolving adaptive networks as 
well. We define the dynamical rules of our pair-based proxy model to follow those of the original network 
model as closely as possible, based on the typical changes that occur in the network model during infection, 
recovery, birth or rewiring events. 

Our approach relies on an approximation that neglects long-range correlations within the network, equivalent 
to the usual PA which is often employed when writing down deterministic effective descriptions. As such, 
the behaviour of the PBP may only be considered an accurate representative of the full adaptive-network 
model in situations where typical states of the network are not too spatially heterogeneous. The agreement 
is aided in our model by the presence of a low-level random rewiring process which acts to 'stir' the network. 
As this rate is increased we observe a corresponding improvement in the accuracy of the predictions made 
by the PBP. As well as improving the predictive power of our methods, we would argue that the presence of 
background of random rewiring is not unrealistic as we expect that most naturally occurring networks are 
subject to extrinsic factors (e.g. movement of agents) which randomise their structure. 
Comparisons with network simulations demonstrate that the PBP is able to describe the behaviour of the 
underlying adaptive-network model on the aggregate deterministic level (fixed points are reproduced to a 
good accuracy), and crucially also on the level of Gaussian fluctuations for large, but finite networks. To 
our knowledge no previous (semi-) analytical work exists which would predict the power spectra of quasi- 
cycles in adaptive-network models to the accuracy shown in Fig. [TJ While our work focuses on the specific 
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application to models of epidemic spread we expect that they can be transferred to adaptive-network models 
of other processes as well, in particular for example models in game theory [2] [19j [20] , opinion dynamics [21] 
or ecology [T]. 

Almost as a by-product we also obtain access to oscillatory behaviour of the network itself. In the adaptive- 
network model the dynamics on and of the network co-evolve, and as such one would expect the network 
structure itself to oscillate in regimes in which the disease dynamics shows quasi-cycles. This is indeed 
confirmed by simulations, in which we observe noise-driven oscillations with amplitude proportional to 
1/y/N. Moreover, our analysis using the PBP allows us to obtain accurate results for the spectral properties 
of these oscillations in network structure. Since our approach is analytical, we are able to systematically 
explore parameter space, for example to determine the relationship between the character of the oscillations 
and the rates of infection and smart rewiring. We expect that other adaptive-network models will also exhibit 
fluctuations in their structure which can be captured within our approach. 
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